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Abstract: We discuss the question of the relevance of perturbative QCD calculations 
for analyzing the properties of the dense medium produced in heavy ion collisions. Up 
to now leading order perturbative estimates have been worked out and confronted with 
data for quenched large p±_ hadron spectra. Some of them are giving paradoxical results, 
contradicting the perturbative framework and leading to speculations such as the formation 
of a strongly interacting quark-gluon plasma. Trying to bypass some drawbacks of these 
leading order analysis and without performing detailed numerical investigations, we collect 
evidence in favour of a consistent description of quenching and of the characteristics of the 
produced medium within the pQCD framework. 
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1. Introduction 

There is convincing evidence that hadron spectra at high p± in nuclear collisions are 
strongly suppressed [I], |2|. For neutral pion production in Au-Au collisions at RHIC the 
nuclear modification factor Raa(p±) is measured to be about 0.2 — 0.3 in the range of 
transverse momenta as large as 10 GeV/c < p± < 20 GeV/c [Q, |J, |5[. This quenching is 
usually attributed to radiative parton energy loss, for which the relevant expressions are 
obtained from perturbative QCD calculations || fjj. They are used to extract and ana- 
lyze the properties of the deconfined medium produced in the collisions. As a consistency 
requirement these properties should be compatible with the perturbative framework. 

However, the detailed analysis given in || [| results in a (time-averaged) value for 
the transport coefficient q which characterizes the medium, exceeding 5 GeV 2 / fm, almost 
a factor of 10 (or even more) bigger than a typical perturbative estimate at the energy 
density expected for ^snn = 200 GeV Au-Au collisions! 

An independent work [1C], based on [11], calculating the quenching factor for hard 



pion production and comparing it with data measured at RHIC obeys this requirement of 
perturbative consistency: an averaged value of q ~ 0.3 — 0.4 GeV 2 / fm is found, corre- 
sponding to an energy density e ~ 2 GeV/fm 3 , as expected from pQCD for a deconfined 
equilibrated plasma. However, this value is based on imposing arbitrarily a mean path 
length for the jet of about L = 5 fm, whereas for the denser medium discussed in || a 
characteristic path length of L ~ 2 fm is obtained. 

A related work on the nuclear modification factor for leading large p± hadrons (pions) 
1 12], assuming a thermalized medium (and also taking into account absorption of thermal 



partons), shows that a value of Raa is obtained which is compatible with measurements 
and the perturbative framework. This work based on the AMY [O] formalism describes 
the coherent gluon radiation, incorporating the Landau-Pomeranchuk-Migdal (LPM) ef- 
fect in the range of gluon energies uj > ujbh, where ujbh (sometimes denoted by Elpm) 
corresponds to the transition energy to the incoherent Bethe-Heitler radiation regime. 

The analysis in ||], following the explicit calculations of "quenching weights" [14], does 
not distinguish between these two regimes. On the other hand, it does impose a kinematical 
constraint, taking properly into account the effect due to the transverse momentum phase 
space of the emitted gluon. This effectively constrains the soft LPM emission by imposing 
a lower energy cut-off oj, which depletes the gluon energy distribution. 

In this note we argue that the introduction of this cut-off uj is a priori not accurate 
enough because, as we shall see below, uj is (much) smaller than ujbh and thus not relevant 
for the LPM regime. This actually is one crucial reason for the large value of q found in 
Q. We instead show that using ujbh as the proper cut-off for the validity of soft LPM 
gluon emission may lead to values of q compatible with perturbative estimates. 



The line of our arguments treating radiative energy loss follows the BDMPS [15, 16] 
- Zakharov [17, [l^] - Wiedemann [p~9| ] approach, as it is reviewed in || 0]. The trigger 
bias induced by the steeply falling large p± vacuum production cross section of produced 
hadrons/neutral pions is treated as described in pOf| . 

The infrared sensitivity of the quenching factor has already been commented upon in 
[20], where it is emphasized that the energy ujbh plays a central role. 

2. Limits on the quenching factor 



As it is discussed in rather great detail in [|20| the quenching effect is expressed by the 
factor 
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where it is justified to express the probability D(e) for emitting the energy e into the 
medium by a Poissonian energy distribution 
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This expression is assumed to be valid for the emission of soft primary gluons. In the LPM 
regime uj > ujbh the bremsstrahlung spectrum dl /duo is given in Correspondingly the 
multiplicity of LPM gluons with energies larger than u is given by 
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For uj significantly less than the characteristic energy [ITj 
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but larger than to bh > the number of gluons is well approximated by 
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The following remarks are crucial for the subsequent analysis: 

• The probability D(e) is normalized by J deD{e) = 1. 

• The ratio of cross sections in ( |2.1[) is well approximated by 

da™ cuvm (p ± + e)/dp 2 l ( p± \ n ( ne\ 



exp , (2.6) 
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when expressed in terms of an effective exponent n, which may depend on n 
n(p±). In the following the approximation (|2.6|) is used. 



• Concerning the underlying parton interactions one has to distinguish quark versus 
gluon jets. Since our concern is mainly neutral pion production at RHIC in the 
range 10 < p± < 20 GeV/c, we effectively assume a dominating quark jet, radiating 
off (soft) gluons, therefore we take Cr = Cf = 4/3 in (|2.5| ). This assumption is 
supported by the analysis in ||. 

• The quenching factor Q(p±) corresponds to the experimentally measured ratio || 

RAA ^ = 7 M > W 

< N co u > dN NN 

for central A- A collisions versus nucleon-nucleon (NN) collisions. In the following 
only neutral pion production at pseudo-rapidity rj = is considered. 

• The transverse momenta of the produced pions are not asymptotically large, although 
only leading order pQCD calculations are considered. 

Neglecting in the following the contribution from Bethe-Heitler emission (Appendix 
A), the quenching factor due to LPM emission becomes 

Q(p±)^J^ deD{e) expj-^j = expj-^ ^-l-exp(-^ dwj, (2.8) 

where we take as a lower cut-off the energy lubh, which indicates the transition between the 
Bethe-Heitler and the LPM regime. By partial integration and with ( |2.3[ ), the quenching 
factor becomes 

Q(p±) = exp < —N(lobh) 1 - exp ( - nuJBH | l-expi — — / N(uj) exp ( -— ) dtu \ 
I L V P± )\) I P± Jljbh V P±J J 

(2.9) 

As N(u) decreases with increasing gluon energy pQ ], one finds the following two bounds 
for Q{p±): 

Qmin(p±) = exp [-N{u BH )} , (2.10) 



i.e. by replacing N(u) = N{ubh) in the integrand of the integral in (2.9), and 



Qmax(p±) = exp < -N(u; B h) 
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i.e. by neglecting the second exponential factor of fl2.9| ). So that, the experimental ratio is 
required to be bounded as follows, 

Qmin (p±) < RAA(P±) < Qmax(P±) ■ (2.12) 

One notes that actually Qmi n (p±) does not depend on and that Qmax (p±) approaches 
1 for large p± » ulobh- For a fixed value of p±/n ~ 0(1 GeV), the interval Qmax — Qmin 
becomes rather tight, when u>bh is not a small energy. 

When comparing the quantities Q m in,max with the ones defined in [||, |14|] Qmin is 
related to the discrete weight of the probability distribution D(e) by po = Qmin-, where, 
instead of the soft cut-off Cj the Bethe-Heitler one ojbh has to be taken, as required by 
the LPM radiation formalism. According to our analysis Qmax is the relevant, important 
quantity, related to the continuous part of D(e), which translates to Qmax — Qmin- 

One may notice that the suppression factor Q(p±) of ( |2,9D does satisfy (for the same 
value of uj c ) 

Q(p±,ubh)>Q(p±,Cj), (2.13) 

where the value of the lower cut-off is indicated. But here is where Qmax comes into the 
game. Its actual value and the ones of the various energy scales cannot be left out of the 
discussion ! 

Typical estimates of the medium characteristics as constrained by the experimental 
results are discussed in the next section. In particular we find that ojbh ~ 1.5—2.0 GeV. As 
a consequence, for values of p±/n ~ 1 GeV, the bounds Qmin, Qmax are indeed constraining. 
On the contrary, taking Cj as the lower energy cut-off, one finds that due to Cj « ujbh, 



Qmin and Q m ax differ significantly and the relation ( 2.12 ) is no longer constraining. 

Replacing ujbh by Cj in the estimate for Qmax, the constraints are not very tight indeed: 
as an example when choosing Cj ~ 0.45 GeV, guided by relation (3.6) to be derived in the 
next section, we find Qmax — 0.55, instead of Q max — 0.30. In both cases Q m in = 0.2. 



3. Kinematical and consistency constraints 

As mentioned earlier, the detailed discussion given in g] in order to determine the medium 
induced gluon radiation intensity dl/duj takes into account the kinematical constraint as- 
sociated to the transverse momentum phase space of the emitted gluon. This constraint is 
not implemented in earlier works [15]. The effect of the kinematical limitation is obtained 
by estimating the ratio k±/uj in the LPM regime: in this coherent regime the transverse 
momentum k± of the emitted gluon may be given by 



where /i is the typical transverse momentum transfer in a single scattering (i.e. the Debye 
mass screening the gluon exchange) and X g the gluon mean free path, such that t co h/X g 
is the number of coherent scattering centers in the medium which the gluon encounters 
before being emitted after the time t m h — 1ujjk\. One finds 



k\ ~ , (3.2) 

where q ~ [i 2 /\ g . As a consequence, since k± < u, gluons have to be emitted dominantly 
above the energy u) defined by 

u = (2qf/ 3 =u J J-) , (3.3) 



R) 

where it is convenient to introduce the dimensionless parameter [^l], |2^, |l4| 

R = u; c L= ^qL 3 . (3.4) 

Now, we should take into account that the multiple scattering formalism used through- 
out the derivation of dl/dto requires the condition [23, 15 1 



Xgfi » 1 . (3.5) 

Using the fact that lubh may be expressed as ujbh — ^g^ 2 we obtain the following para- 
metric inequality 

a 2 1 ' 3 
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As a consequence of this inequality, the Bethe-Heitler energy remains the proper and 
relevant lower limit for the validity of the LPM gluon emission spectrum. 

Let us summarize what is already known about the implementation of the various 



above mentioned constraints. In the BDMPS framework fig , 16] where no kinematical 
cut-off is imposed on the k± integration, i.e. in the limit R — > do, the gluon number only 
depends on the ratio uj/uj c : N(uj) = N{uj/u) c ). As a consequence the resulting quenching 
factor Q(p±) is effectively a scaling function in the variable X = p±/(nu> c ) (20], such that 
in the relevant analysis given in [ |10| only the characteristic gluon energy uj c is extracted 
from the the neutral pion single-inclusive data measured by the PHENIX Collaboration 
in Au-Au collisions |3|, 0], and found to be u c = 20 — 25 GeV. Moreover, as already 
mentioned, the path length is arbitrarily chosen so that the medium characteristics are 
not quantitatively constrained. In || the typical value of the parameter R relevant for 
the description of RHIC data on pion production at ^snn = 200 GeV is estimated to be 
R ~ 1000, equivalent tow~ 0.016 oj c . 



4. Comparison with experiment 

We shall first discuss a few semi-quantitative estimates of the parameters describing the 
medium extracted from the comparison with data within the framework described in the 
previous sections. 



4.1 Averages 

We concentrate on the data for Raa(p±) with p± > 10 GeV/c: Raa(p±) — 0.2, remaining 
essentially flat up to p± ~ 20 GeV / c. 

In order to start the discussion, we use these data to obtain values for Q m in and Qmax, 
corresponding to the experimental error bars. Taking Q m in — 0.2 leads to N{ujbh /w c ) — 
1.6, which allows to estimate ujbh/uj c — 3.5 • 10 -2 when using ( ^ ) and a s = 1/2 (see also 
Fig. 1). When we take p±/n ~ 0(1 GeV), observing that the experimental error bars 
allow us to fix the value of Q m ax < 0.3, we deduce the typical value of ujbh- From ( 2.11| ) 



we obtain ujbh > 1-4 GeV, in agreement with thermal estimates (Appendix B). Taking 
ubh — 1-6 GeV , we find u> c ~ 45 GeV. We note that this value is bigger by a factor 2 
than the one extracted in [|i~0| ]. 

Depending on the typical average path length we derive estimates for the time — averaged 
transport coefficient, 

„ _ 2uc _ 18 GeV 2 

q ~U~ {L[fm}) 2 ~Jm~ ' ( ' J 

i.e. q ~ 1.1 GeV 2 /fm for L = 4 fm, and g ~ 2.0 GeV 2 / fm for L = 3 /m. This 
corresponds to values of R ~ 900 — 700, not far from i? ~ 1000 in |8|. 

We observe that these estimates are indeed valid beyond the R = oo limit, since in the 
region of ujbh/uj c > 3.5 10~ 2 the sensitivity on R, for R > 1000, becomes weak. 

The cut-off ujbh is effectively of order 1 — 2 GeV whereas u>, as imposed by Eskola et al. 
[8] is also 1 GeV, when using relations given in [14]. One may wonder then why choosing 
one or the other is a crucial feature, independently of judging the validity of using the LPM 
spectrum away from to > ujbh- One way to understand qualitatively this fact is that in 
the analysis of Ref. [14] , the u cut-off is implemented effectively in the k± integration for 
the emitted gluon whereas here the ujbh cut-off appears as an IR cut-off on the emitted 
gluon energy in the expression of Q(p±). That has consequences on the effectiveness of the 
cut-off. 

The question of the impact on the value of q of choosing uj - instead of ujbh - as the 
cut-off for the spectrum, should now be discussed. Indeed the above indicates that although 
larger than expected from leading order estimates based on the presence of a thermalized 
and ideal QGP (Appendix B), the values obtained above for q are much smaller than the 
ones quoted in (8|, namely 5 < q < 15 GeV 2 / fm. 

An even smaller value may be obtained with L = 5 fm, namely q ~ 0.7 GeV 2 / fm. In a 
thermal gluonic system this implies an (time — averaged) energy density of e ~ 4 GeV/ fm? . 

It is obvious that the values of q and L are strongly correlated, namely a large transport 
coefficient, corresponding to a dense medium implies a shorter path length L, and vice 
versa. For a realistic average path length of L ~ 3 fm in the case of Au-Au collisions under 
consideration the prefered value of the time — averaged transport coefficient becomes 

q ~ 2 GeV 2 /fm , (4.2) 



which may still be accomodated into the pQCD framework, at least within the uncertainties 
of LO approximations, contrary to the "strong" QGP values of ||]. 
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Figure 1: The gluon multiplicity as a function of u>, for R = oo (solid curve) and R — 1000 (dashed 
curve) taken from E3], but for a s = 1/2. 



In fact, the actual values of the cut-off cannot be left out of the discussion. Taking 
numbers quoted above: q = 10 GeV 2 /fm and L = 2 fm, leads to u) c = 100 GeV. If as 
indicated above, we take uj/uj c — 10~ 2 , one finds cj ~ 1 GeV, and thus correspondingly 
ujbh — 3 — 4 GeV, which is too large to make sense for energies/transverse momenta under 
consideration! If, on the other hand, we want to have a reasonable value oIuibh — 1-4 GeV , 
from the start, imposing correspondingly that the value of u is a factor 3 — 4 smaller and 
keeping uj/uj c ~ 10 -2 , we find ui c and thus q, 3 -4 times smaller. We use the value of 



Qmax < 0.3 and Q 



0.2 as a way to constrain the relevant parameters: from Q max , we 



take the cut-off to be ~ 1.4 GeV. This cut-off can only be ojbh- Then from Q m in — 0.2, 
we deduce u c . Finally, we find, depending on the length L, reasonably small values of q. 
A determination of the average L should be possible with the expression for Q m i n (p±), 



(|2.10[) , which depends on L. Since ujbh <C ui c , we use to a good approximation ( [2.5D to 
determine N(lobh) = N(ubh/w c ) and insert (neglecting logarithmic factors) 
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defining the quark mean free path A„ 



C F n 9 



\ q , such that 



(4.3) 



<3min(P±) ^ exp 



2a s C F 



IT 



Cf \ 



+ In 2 In 



L 2 iV c 2 



1.44 



(4.4) 



In leading order L/X q S> 1 the dependence with respect to the path length has the 
typical characteristic behaviour of a survival probability exp [— L/A e //], where A e // — 
2a N \ - \ for a s ~ 1/2. Without further geometrical restrictions the mean path length 
(L) would be just given by the mean free path of a quark jet in the medium, (L) ~ X q . 



A better estimate of L/X q , however, is obtained by taking the full expression ( |4.4[ ) into 
account: for Q m in — 0.2 the corresponding ratio is L/X q ~ 3.35. 

In order to obtain L, we estimate the mean free path X q = 9/4A 9 from q and u>bh by 



with q ~ 2 GeV 2 / fm and ojbh — 1-6 GeV and find L ~ 3 /m, consistent with the prefered 
value given above. Note the related estimate for the mass fi: fj, ~ 0.9 GeV. 

4.2 Nuclear geometry 

So far we have considered averaged values for Qmin, Qmax without taking the nuclear 
geometry explicitely into account. We now present a more detailed discussion which, as 
we shall show, leads on a firmer basis to similar conclusions as above as far as the medium 
parameters are concerned. We assume head-on collisions and essentially cylinder-like Au 
nuclei. The quark jet is produced in Au-Au collisions at mid-rapidity and propagates in 
the transverse plane. Following H, one starts from the geometrical transverse path length 



L g eom{s) = -s cos c/) L s + J s 2 cos 2 (j) LS + R 2 Au - s 2 , (4.6) 

where the position at which the parton is produced is denoted by the vector s in the trans- 
verse plane, (pLS is the angle of propagation with respect to this vector. This geometrical 
picture allows us to obtain an average value for {Qmin) by calculating 

in (A/,, " _ Jd 2 sexp{-N(L scom )} 

nK Au 

with |s| < Rau-, the radius of the Au nucleus. In order to obtain N(L geom ) we use (^Jj) 
with uic/uiBH = q/(2uj B H) L 2 geom . 

In Fig. H we plot {Q m in) as a function of (j/lubh, and observe that (Qmin) ~ 0.2 for 

A ; , . GeV (A o , 

q/uBH - 1-4 ■ (4.8) 

fm 

In the same Fig. |2] we plot (Qmax), obtained analogously to (|4.7| ), for different values of 
ojbh (and for n/p± = 1 /GeV). 

In order to have {Qmax) < 0.3 - together with (Q m in) — 0.2 - we find 0.75 < lobh < 
2.0 GeV. 

In Fig. pi we show (Qmin), but this time as a function of A g , obtained from 



{Qmm[Xq)) = , (4.9) 



J d S Qmin(Lg eora / ' \q) 

after inserting L = L geom in (|4.4| ). This way we find X q ~ 0.85 fm (X g ~ 0.38 fm), when 



Qmin — 0.2, in good agreement with (4.5). 

From Fig. ^ we deduce that 0.65 < /i < 1.1 GeV. Finally, within these bounds the 
transport coefficient q ~ A t2 /A g becomes 

q ~ 1.0 - 3.0 GeV 2 /fm . (4.10) 
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Figure 2: (Q m in) (solid curve) and (Q ma x}, respectively, as functions of q/u>BH according to (4.7) 
(Qmax) for different values of lubh :0.75 (dotted), 1.15 (dashed-dotted) and 2.75 GeV (dashed curve) 



Averaging with respect to the nuclear geometry in this straightforward manner leads to a 
final estimate which is compatible with ( [4.2|) . The "uncertainty" given by ( 4. 10| ) may be 
considered as "theoretical error" on the derived medium parameters. 

Concerning the robustness of the estimate (4.10) for q, one has to keep in mind that 
besides the IR sensitivity under discussion, all the estimates are based on LO QCD calcu- 
lations. Nevertheless, the explanation of quenching as being due to the LPM radiation in 
a perturbative regime appears to be a robust statement. 
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Figure 3: {Q m in) (solid curve) and (Q ma x), respectively, as functions of A g according to fl4.S| ), 
(Qmax) for different values of the screening mass fi: 0.65 (dotted), 0.8 (dashed-dotted) and 1.1 GeV 
(dashed curve). 



A value of q ~ 1.8 GeV 2 / fm would correspond to a temperature of T ~ 375 MeV, 
corresponding to an average energy density of e ~ 12.5 GeV/ fm? , in agreement with the 



results quoted in [12|. 

In order to obtain the actual value of q at the very early stage of the collisions, at times 
r ~ l/pj_ < 2 • 10~ 2 /m, one has to include the effects due to the longitudinal [24, 22], but 
also transverse [^] expansion of the dense system during the time of O(L), the jet takes 
to propagate through this medium. 
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A. Estimate of the Bethe-Heitler and absorption contributions 

In the plausible context, where the medium is thermalized, let us investigate how the 
contributions of the Bethe-Heitler and the absorption spectra modify the analysis discussed 
so far in the paper. 

In LO pQCD the Bethe-Heitler-Gunion-Bertsch [26] differential spectrum for inclusive 
gluon production in a medium of length L in the presence of L/X g scatterers reads 
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which after pj_-integration becomes 
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(A.l) 



(A.2) 



Because of the presence of the (non-perturbative) IR-cut, w cu t, this BH-intensity is not 
precisely determined. Nevertheless, an estimate of the gluon energy ujbh may be obtained 
from matching ( |A.2j ) with the LPM expression for the intensity at uj = ujbh, 
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valid for uj -C lu c . When logarithmic factors are taken to be of O(l), i.e. ujbh — 1-65 w cu t, 
we find 

uj BH ^ /U 2 A 9 , (A.4) 

in terms of the screening mass \x and the mean free path of the gluon X g . 
The resulting quenching factor becomes 
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and expanding the integrand in the limit of small u>, we find with ( |A.2j ) 
Qbh(p±) > exp 



a s Cp ( L\ n s 
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Following fl^l , we consider the absorption contribution in the presence of a heat bath. 
It is assumed that for the radiation energy 00 < the spectrum is approximated by 
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i.e. the Bethe-Heitler spectrum multiplied by a Boltzmann factor with temperature T. 
The corresponding quenching factor Q a bs{PL) then becomes 
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The integral may be approximated by 



leading to 
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Finally, multiplying the two quenching factors Qbh(pl) and QabsiPL), leads to a lower 
bound 



QBH{PL)Qabs(PL) > exp <^ + 
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Taking as typical numbers: ujbh — 1-6 GeV, uo cut ~ 1 GeV, T = 350 MeV, p± 
10 GeV, n = 10, and L = 3 fm, X g = 0.38 fm, we find 



QBH{PL)Qabs{PL) > 0.7 . 



(A.12) 



This indicates that we may as a first guess, as already suggested in [12], neglect alto- 
gether the Bethe-Heitler and absorption processes, meaning that the values we thus obtain 
for q are in fact upper bounds, and therefore comforting the perturbative framework. 



B. Thermal parameters 

Let us consider an equilibrated system (N c = 3, Nf = 0) in the weak coupling QCD limit 
1 27] and give a short summary of the elements which enter our analysis. 
In LO given the temperature T the screening mass is 



H 2 = 4vra s T 2 . 



(B.l) 



The gluon mean free path X„ is expressed in terms of the gluon density 
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2C(3)T 3 , C(3) = 1.202 , (B.2) 



and the (transport) gluon-gluon cross section (to logarithmic accuracy) 



when neglecting logarithmic dependence, i.e. 



ln(l/a s )~ — (B.3) 
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1/A fl = ~ -^((3)a s T ~ 2.2 a s T . (B.4) 

7T 

The mean free path for a quark is \ q = 9/4A 9 . The corresponding energy density of this 
gluonic system is 

e = ^T\ (B.5) 
15 V ; 

Typical orders of magnitude may be given e.g for a temperature of T = 400 MeV 

and a coupling a s = 1/2: the screening mass is // ~ 1 GeV, the mean free path A g ~ 

0.45 fm (X q ~ 1 /m), implying ujbh — — 2.25 GeV . The transport coefficient is 

estimated as q ~ /i /A p ~ 2.2 GeV 2 / fm, leading to an energy density of e ~ 17 GeV j fm 

when using q ~ 2 e 3//4 p8| . 
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